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(57) Abstract: An apparatus 10 and method 
3S provided that includes a spatial array of 
unsteady pressure sensors 15 - 18 placed at 
predetermined axial locations Xi - x>r disposed 
axiallyalong a pipe 14 for measuring at least 
one parameter of a solid particle/flnid mixture 
12 flowing in the pipe 14. The pressure sensors 
15-18 provide acoustic pressure signals Pi(t) 
- PN(t) to a signal processing unit 30 which 
determines the speed of sound an,i3t(a)) of 
the particle/fluid mixture 12 in the pipe 14 
using acoustic spatial array signal processing 
techniques. The primary parameters to be 
measured include fluid/particle concentration, 
fluid/particle mixture volumetric flow, and 
particle size. Frequency based sound speed 
is detennined utilMng a dispersion model 
to determine the parameters of interest the 
calculating the at least one parameter uses an 
acoustic pressure to calculate 



wo 03/091671 Al liliillliliilllililiiiililll 



European patent (AT, BE, BG. CH, CY, CZ, DE, DK, EE, 
ES, H, FR, GB, GR, HU, IE, IT, LU, MC, NL, FT, RO, 
SB, SI, SK, TR), OAPI patent (BE, BJ, CP, CG, CL, CM, 
G A. GN, GQ, GW. ML. MR. NE. SN, ID, TG). 

Published: 

— -with international searck report 



— before the expiration of the time limit for amending the 
claims and to be republished in the event of receipt of 
amendments 

For two'letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



wo 03/091671 



1 



PCT/US03/12956 



Apparatus and Method for Measuring Parameters of a Mixtiire Having Solid Particles 

Suspended in a Fluid Flomng in a Pipe 

Cross-Reference To Related Patent Applications 
5 TMs application is a continuation in part of U.S. Patent Application No. 10/349,716 

(Cidra Docket No. CC-0579) filed on January 23, 2003, wMch claimed the benejBt of U.S. 

ProvisiondApplicationNo. 60/351,232 (CidxaDocketNo. 

2002^US. Provisional AppHcation No. 60/359,785 (Cidra 

February 26,2002; U.S. Provisional AppUcationN^^ 
10 0468), filed April24, 2002; U.S. Pro\dsional AppUcation N^o. 60/425,4^ (CidraDocket 

No. CC«-0538), ffledNoyember 12, 2002; and U.S. Prpvisioiial ApjpHcationNo. 60/426,724 
: (Cidra Docket Np. CC-0554), filed November 15, 2002; and is a continuationinpart of U.S. 

Patent Application No. 10/376,427 (Cidra Docket No. CC-0596) filed on February 26, 

2003, wMch claiined the benefit of U.S. Prpvisiond^ 
15 Docket No. CC-0403), filed February 26,2002; and is a continuation-in-part of U.S. Patent 

AppUcation Serial No. 10/349,716 (Cidra Docket No. CC-0579), filed Jan^^ 

which claims the benefit of U.S. Provisional Application No. 60/351,232 (Cidra DocketNo. 

CC-0410), filed January 23, 2002; U.S. Proyisi^^ 

Docket No. CC-0403), filed Febniary 26,2002; U.S. Provisional Application No. 
20 60/375,847 (Cidra Docket No. CC-0468), filed April 24, 2002; U.S. Provisional Application 
No. 60/425,436 (Cidra Docket No. CC-P538), filed Novem^^ 

Provisional Application No. 60/426,724 (Cidra Docket No. CC-0554), filed November 15, 
2002„ all of which are incorporated herein by reference in their entirety. 

25 Teclmieal Field 

This invention relates to an apparatus for measuring the flow passing within a pipe, 
and more particularly to an apparatus and method for measuring the speed of sound 
propagating in the flow, having particles suspended within a continuous fluid, to determine 
parameters, such as particle/fluid ratio, particle size and volumetric flow rate of the flow in 
30 pipes usiag acoustic dynamic pressures. 
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Background Art 

This inventioii provides a method to measure parameters of a fluid/particle mixture 
in a pipe that can be used in many applications, such as in chemical, pharmaceutical, 
petrolemn and power generation industries. In particular, the invention provides a method 

5 to measure pulverized coal and air mixtures used in pulverized fuel delivery systems in 
place in a large percentage of cod. fired boilers used in the power generation iudustiy. 

Currently, well over 50% of the electricity ia the US is generated with coal. While 
coal is considered a cost effective, abundant resource in the US, the use of coal has been 
restricted due ia large part to environmental concerns. To mitigate this impact, the US 

10 Department of Energy and the Power Generation industry have large prograins desi^^ to 
develop technology to reduce the environment effects of burning coal. These Clean Coal 
Initiatives include technology designed to develop improvements in the combustion process 
to improve efficiency while reducing pollutants such as unbumed carbon, ash, and nitrous 
oxide CNOx). 

15 The ability to measure the flow rate and composition of the air /coal mixture withm 

the cp4 pipes is an important aspect of any sy stem or strategy 

performance of the PF delivery system. The industry recognizes ttds and therefore has been 
developing a wide variety of technologies to perform this mea^ These include 

probe based and sampling devices, as well as real time meters based on a wide variety of 
20 technologies including electrostatic charges, microwaves, and ullxasonic. 

Summary of the Invention 

Objects of the present invention include providing a system for measuring the speed 
of sound propagating through a particle/fluid mixture in pipes in iudustrial boiler systems 
25 and related processes, such as coal fned boiler systems, to determine particular parameters 
of the mixture. 

According to the present invention, an apparatus for measuring at least one 
parameter of a particle/fluid mixture in a pipe iacludes a spatial array of at least two 
pressure sensors, disposed at different axial locations along the pipe. Each of the pressure 
30 sensors measures an unsteady pressure within the pipe at a corresponding axial location. 
Each of said sensors provides a pressure signal indicative of the unsteady pressure within 
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the pipe at said axial location of a corresponding one of said sensors. A signal processor, 
responsive to said pressure signals, provides a signal indicative of the at least one parameter 
of the mixture in the pipe. 

According to the present invention, a metiiod for measuring at least one parameter of 

5 a particle/fluid mixture in a pipe includes measuring unsteady pressures within the pipe at at 
least two predetemuned axial measurement locations along the pipe to provide a pressure 
signal indicative of the unsteady pressure within the pipe at each of the at least two 
predetenrdned axial measurement locations. Further the method iacludes calculating the at 
least one parameter of the particle/fluid mixture in the pipe usiag the unsteady pressure 

10 measured at the axial measurement locations. 

The foregomg and other objects, features and advantages of the present invention 
will become more apparent in Hght of the following detailed description of exemplary 
embodiments thereof. 

» 

15 Brief Description of the Drawings 

Fig. 1 is a block diagram of a flow meter for measuring the speed of sound of the 
fluid/particle mixture jaowing with apipe, in accordant 

Fig. 2 is a schematic diagram of a pulverized fuel (PF)/air mixture parameter 
measurement system within a coal fired boiler system, in accordance with the present 
20 invention. 

Fig. 3 is a magnified photograph showing particle size of coal typical of the system 
shown in Fig. 2. 

Fig. 4 is a plot of the speed of sound of a mixture versus the firequency in air/coal 
mass flow ratio, in accordance with the present invention. 
25 Fig. 5 is a plot of actual data and a model of the speed of sound as a function of 

frequency for air/coal mixtures, in accordance with the present invention. 

Fig. 6 is a plot showmg the standard deviation of sound speed versus frequency for 
various arrays of a PF/air mixture parameter measurement system, in accordaace with the 
present iavention. 
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Fig. 7 is a plot of sound speed as a ftmction of firequency for air/coal mixtures with 
fixed particle size (50 mm) and varying air-to-fael mass Ratio in accordance with the 
present invention. 

Fig. 8 is a plot of sonnd speed as a function of firequency for air/coal nnxtures with 
varying particle size where the air-to-fiiel mass ratio is equal to 1.8 in accordance with the 
present invention. 

Fig. 9 is a plot of sound speed as function of ak/coal ratio in accordance with the 
present invention. 

Fig. 10 is a flow diagram ofan optimization procedure employed to determine a^^ 
to-fuel ratio and particle size from analytical model and experimentally detemiined 
dispersive speed of sound data in accordance with the present invention. 

Fig. 1 1 is a plot of the results of the optimization procedure of Fig. 10 applied to 
data recorded from an array of sensors listening to flow in a six inch circular duct, 50 pm 
particle size, 100 ft/sec air flow rate with an air-to-fuel ratio of 1 

Fig. 12 is aplot of the results of the optimization procedure of Fig. 10 applied to a 
series of data sets with varying air-tp-fiiel ratio. 

Fig. 13 is ako) plot of data processed from an array of pressure sensors use to 
measure the speed of sound of a coal/air anixture flowing in a pipe, in accordance with the 

present invention. , 

Fig. 14 is a side elevational view of a plurality of pressure sensors, having PVDF, 

clamped to the outer surface of tihe pipe, in accordance with the present invention. 

Fig. 15 is a partial perspective view of one of the pressure sensors of Fig. 14. 

Best Mode for Carrying Out the Invention 

Referring to Fig. 1, a flow meter 10 embodying the present invention is provided 
that measures a number of parameters/characteristics of a mixture 12 of solid particles 
suspended within a continuous fluid flowiag withm a pipe or conduit 14, whereia a fluid is 
defined as a liquid and/or a gas. The flow meter may be configured and programmed to 
measure the speed of sound propagating through the mixture. The flow meter can measure 
at least one of the following parameters of the nrixture flow 12: the fluid/particle 
concentration (volumetric phase fraction), the volumetric flow rate, the size of the soUd 
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particles, .the mass flow of the m&ture and the velocity of the mixture. To determine any 
one of these parameters, the flow meter 10 measures the unsteady pressures created by the 
speed of sound (SOS) propagating through the mixture flowiag in the pipe 14, which will be 
described in greater detail hereinafter. 

5 The solid particles of the mixture 12 inay be of any size, shape and materiaL For 

example, the particles may be small in size as in the form of a po-v^der, in a graaular form, 
or greater in size. The flow meter 10 can be used in any appUcation that carries solid 
particles suspended in a fluid through a pipe, such as in chemical, pharmaceutical, 
petroleum aad power generation appHcations. For example, the present invention is well 

1 0 suited to measure the parameters (e.g. ahr/coal ra.tip, particle size) for power generation 
systems that use pulverized coal to fire the jRimace a steam boiler system. 

As one example, the present invention will be discussed in the context of a 
Pulverized Fuel (PF) delivery system for power generation, but one wiU appreciate that the 
flow meter can be applied to any number of other applications, as discussed hereinbefore. 

15 A representative PF delivery system 1 is shown in a coal fired boiler system 2 in Fig. 2. 

• ■ . ^ ...... ... 

The coal is pulverized in a iniU 3 and entrained m 

fan4 to transport the PF/air mixture ^ 12 for dieUyery to the Typical 

furnaces can have >50 coal pipes, each 12-20 inch in diameter. T?ypically, a large utility 
boiler >300 Mw, can have 4-1 1 pulv^ing miUs feeding the furnace. The ability of the PF 
20 dehvery system to deliver the proper amount of fuel and air to the furnace through these 
multiple cod pipes, both coUectively and individuaUy, has a strong M 

performance and emissions jfrom the coal fired boiler. 

As is known, non-uniformities in the PF delivery system 1 can result in variation of 
the fuel to air ratios, causing hot spots, regions of high NOx generation, and unbumed fixel. 
25 The connection between performance of a PF fiiel delivery system 1 and a boiler system 2 is 
well recognized. The flow meter 10 embodying the present invention is capable of 
measuring the fixel to air ratio and particle size of the pulverized coal provided to the 
furnace to thereby provide feedback to the operator to provide more efScient combustion of 
the coal. 

30 As described hereinbefore, the flow meterlO, of the present invention may be 

configured and programmed to measure and process the detected unsteady pressures Pi(t) - 
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Pi^t) created by acoustic waves propagating through the xnixture to determine parameters of 
the mixture flow 12. One such flow meter 10 is shown in Fig. 1 that measures the speed of 
sound (SOS) of one-dimensional sound waves propagating through the fluid/particle 
mixture to determine the composition the mixture, namely the liquid/particle ratio of the 

5 mixture. The flow meter is also capable of deterniining the average size of the pa^^^ 

velocity of the mixture, and the volumetric flow rate of the mixture. It is known that sound 
propagates through various mediruns at various speeds in such fields as SONAR and 
RADAR fields. The speed of sound of a mixture within a pipe 14 may be determined using 
a number of known techniques, such as those set forth in U.S. Patent No. 6,354,147, entitled 

10 'TFluid Parameter Measurement in Pipes Using Acoustic Pressures", issued March 12,2002, 
and U.S . Patent A.pplication Serial No. 10/007,749, entitled "Fluid Parameter Measurement 
in Pipes Using Acoustic Pressures", filed November 7, 2001, each of which are 
incorporated herein by reference. The present invention utilizes at least one flow meter 10 
to determine various parameters of the liquid/particle mixture, wherein one of the 

15 parameters is the speed at which sound travels within the mixture pipe system as will be 
more fully described herein below. 

the nnxture 12 is measured by passively listening to the flow with an array of unsteady 
pressure sensors to detetmine the speed at which pne-dimensional compression waves 

20 propagate through the liquid/pacrticle mixture conta^ 

As shown in Fig. 1 , the flo^ meter 10 has an array of at least three acoustic pressure 
sensors 15,16,17, located at three locations xipc2,X3 axially along the pipe 14. One will 
appreciate that the sensor array may include more than three pressm'e sensors as depicted by 
pressure sensor 18 at location xn. The pressure generated by the acoustic waves may be 

25 measured through holes in the pipe 14 ported to external pressure sensors 15 - 1 8 or by 

other techniques discussed hereinafter. The pressure sensors 15 - 18 provide pressure time- 
vaiying signals Pi(t),P2(t),P3(t),PN(t) on lines 20,21,22,23 to a signal processmg unit 30 to 
known Fast Fourier Transform (F¥T) logics 26,27,28, 29, respectively. The FFT logics 26 - 
29 calculate the Fourier transform of the time-based input signals Pi(t) - pN(t) and provide 

30 complex frequency domain (or fiequency based) signals Pi(c0),P2(o)),P3(ci)),PN(Q>) on lines 
32,33,34,35 indicative of the fiequency content of the input signals. Instead of FFT's, any 
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other technique for obtaining the jarequency domain characteristics of the signals Pi(t) - 
PwCt), may be used. For example, the cross-spectral densily and the power spectral density 
may be used to form a frequency domain traasfer functions (or frequency response or ratios) 
discussed hereinafter- 

5 The frequency signals Pi(cd) - Pn(co) are fed to ajnix-Mx Calculation Logic 38 which 

provides a signal to liae 40 indicative of the speed of sound propagating through the 
raixture amu(i«>), wMch is a fm 

signal is provided to map (or equation) logic 42, which converts ^ia>) to a percent 

composition of the PF/ah* mixture and provides a %Comp signal to liae 44 indicative 
10 thereof (as discussed hereinafter). Also, if tiie Mach number Mx(cd) is not negUgible and is 

desired, liie calculation logic 40 may also provide a signal Mx(co) to line 46 indicative of the 

Mach number Mx((jo) which is a fimction of frequency . 

For ckcular ducts or pipes 12 as shown inFig. 1, only plane waves propagate for 

frequencies below the cut-on frequency (ref Acoustics of Ducts and MufQers, MJ. Mtmjal, 
15 John Wiley & Sons, New York, 198^^^ 

For a mixture with a sound speed of 500 m/sec in an 1 8 inch pipe, the cut-off 
20 frequency is approximately 600 Hz. Thus, for this example, only one-dimensional acoustic 
waves propagate below 600 Hz. It is important to note that one-dimensional waves can still 
propagate above this frequency, but higher order modes may or may not be present. 

More specifically, for planar one-dimensional acoustic waves ia a homogenous 
mixture, it is Icnown that the acoustic pressure field P(x,t) at a location x along a pipe, where 
25 the wavelength X of the acoustic waves to be measured is long compared to the diameter d 
of Ihe pipe 12 (i.e., X/d »1), may be expressed as a superposition of a right traveling wave 
and a left traveling wave, as follows: 

P(;c, 0 = [as -'^-^ + Be ^''^'^^ )b''^' Eq. 1 

where A,B are the fi"equency-based complex amplitudes of the ri^t and left traveling 
30 waves, respectively, x is the pressure measurement location along a pipe, © is frequency (in 
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rad/sec, where a)=27uf), and krM are wave numbers for the right and left travelling waves, 
respectively, which are defined as: 

^ and kf s ^-^ — — ^ ^ ^ ^ Eq. 2 



K 



a> 


1 




l-M^{ai) 



where amix(G>) is the speed of sound of the mixture in the pipe, co is ftequency (in rad/sec), 
5 and Mx(co) is the axial Mach number of the flow of the mixture within the pipe, where: 

where Vmix is the axial velocity of the niuxture. For non-homogenous mixtures, the 
axial Mach number represents the average velocity of the mixture and the low firequency 
acoustic field description remains substantially unaltered. 
1 0 The firequency domain representation P(x,ca) of the time-based acoustic pressure 

field P(x,t) Withm a pipe, is the coefficient of the e*"'4erm of Eq. 1, as foUows: 

P{x,(o)--Ae^'^^Be^^^ Eq.4 

Referring to Fig. 1, we have found that using Eq. 4 for P(x,a)) at three axially 
distributed pressure measurement locations xi,X2,X3 along the pipe 12 leads to an equation 

1 5 for araix as a function of the ratio of firequency based pressure measurements, which allows 
the coefficients A,B to be eliminated. For optimal results, A and B are substantially constant 
over the measurement time and substantially no sound (or acoustic energy) is created or 
destroyed in the measurement section. The acoustic excitation enters the test section only 
through the ends of the test section 51 and, thus, the speed of sound within the test section 

20 51 can be measured independent of the acoustic environment outside of the test section. In 
particular, the firequency domain pressure measurements Pi((o),P2(o>),P3(co) at the three 
locations Xi,X2,X3, respectively, along the pipe 12 using Eq. 1 for right and left traveling 
waves are as foUows: 

Pi(©) ^P{x = x,,ai) = Ae^""^' ^Be^''^''' Eq. 5 

25 PzCco) =P(;c - x^,(D) = ^e^^-^^ +5e"**'M2 Eq. 6 

P3(ct>) =p(jc = x^,(d) = Ae^'^-^' '^Be'"'^'''' Eq. 7 
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10 



where, for a given jfrequency, A and B are arbitrary constants describing the acoustic field 
between the sensors 14,16,18. Fonning the ratio of Pi(a))/P2(G>) from Eqns. 6,7, and solving 
for B/A, gives the following expression: 



B 



Eq. 8 



ik{Xi 



where R is defined as the reflection coefficient 

Forming the ratio of Pi(a))/P3(a)) from Bqs. 5 and 7 and solving for zero gives: 



0 



Eq. 9 



where R=B/A is defined by Eq. 8 and kr and M are related to anux as defined by Eq. 2. Eq. 9 
may be solved nnmerically, for example, by defining an "error" or residual term as the 
magnitude of the left side of Eq.9, and iteratiiig to minimize the error term. 



mag 



Error 



Eq- 10 



The data from the array of sensors may be processed in any domain, including the 
frequency/spatial domain, the temporal/spatial domain, the temporal/waye-number domain 
15 or the wave-number/frequency (k-co) domain. As such, any Imown array processing 
technique in any of these or other related domains nmy be used if desired. 

Also, some or all of the functions within the signal processing unit 30 may be 
implemented in software (using a microprocessor or computer) and/or firmware, or may be 
implemented using analog and/or digital hardware, having sufficient m^ory, interfaces, 
20 and capacity to perform the fimctions described herein. 

Acoustic pressure sensors 15-18 sense acoustic pressure signals that, as measured, 
are lower frequency (and longer wavelength) signals than those used for ultrasonic flow 
meters of the prior art, and thus the current invention is more tolerant to iohomogeneities in 
the flow, such as roping and other time and space domain inhomogeneities within the flow, 
25 even where entrenchment or coal "roping" is imhlcely such as following a bend. The term 
"roping" is a termlcnown to those sldUed in this art which represents a form of severe 



wo 03/091671 PCT/US03/12956 

10 

spatial and temporal mal-distribution induced in mixture flows of widely different 
component densities. It is a condition where a large portion of the coal flow is in a band 
running dong one side of pipe 14. 

In addition, the present invention incorporates the compliance of the pipe 14 to 

5 determine the effective speed of sound of Ihepipe/PF/akinixture system. The acoustic 
pressure signals Pi(t) - PwCt) are generated within the PF/air mixture of the pipe 14 by a 
variety of non-discrete sources such as remote machinery, mills, fans 4 (Fig. 2), valves, 
elbows, as weU as the PF/air mixture flow itself. It is this last source, the PF/air mixture 12 
flowing withinthe pipe 14, which is a generic source of acoustic noise that assures a 

10 niinimum level of acoustics for any PF/air niixture piping systems for which the present 
invention takes unique advantage. The flow generated acoustics increase with mean flow 
velocity and the overall noise levels (acoustic^I^^ pf the 

generating mechanism and the damping mechanism. As such, no external discrete noise 
source is required within the present invention and thus may operate using passive listening. 

15 While the flow meter 10 passively listens to the mixture flow 12, the present invention 

contemplates adding an acoustic source to inject a deske acoustic wave into the flow to be 
measured, such as by compressing, vibrating and/or tapping the pipe, to name a few 
examples. ■ 

For certahi types of pressure sensors, e.g., pipe strain sensors, accekrometers, 
20 velocity sensors or displacement sensors, discussed hereinafter, it inay be desirable for the 
pipe 14 to exMbit a certain amount of pipe compUance. 

Alternatively, to minimize any error effects (and the need for the corresponding 
cahbration) caused by pipe compliance, the axial test section 50 of the pipe 14 along wh^e 
the sensors 15 - 18 are located may be made as rigid as possible. To achieve the desired 
25 rigidity, the thicloiess of the wall of the test section 50 may be made to have a 

predetermined thickness, or the test section 50 may be made of a very rigid material, e.g., 
steel, titanium, Kevlar®, ceramic, or other material with a high modulus. 

The length of the array (aperture) AX of the pressure sensors (15-18) is at least a 
significant fraction of the measured wavelength of the acoustic waves being measured. As 
30 wiU be described in greater detail, the acoustic wavelength to be measured is a function of 
at least the dispersion characteristics of the mixture 12, wherein the dispersion characteristic 
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is a function of at least the size and mass of the particles, and the viscosity of the fluid. The 
greater the dispersion of the mixture (e.g. the greater the size and mass, and/or the less 
viscous the fluid), the longer the length of the array is needed. Conversely, the lesser the 
dispersion of the mixture (e.g. the lesser the size and mass, and/or the more viscous the 

5 fluid), the shorter the length of the array is needed. 

Further, it is within the scope of the present that the spacing of the pressure sensors 
may be known or arbitrary, provided the location of the sensors is known. The sensors IS- 
IS may also be equi-spaced (as shown in Fig. 1) or any noii-even or npn equi-spaced 
location, as will be described in greater detail hereinafter. One will appreciate that as few as 

10. two sensors are required if certain information is Imown about the acoustic prpperti^^ of the 

PF/air mixture piping system. 

As discussed, the flow meter 10 measures the speed of sound of one-dimensional 
sound waves propagating through the fluid/particle mixture to detexmine the composition of 
the mixture. Specifically, the speed of sound propagating through dilute solid/air mixtures 

15 can be directly related to the mass firactionpailicles of the jHow. A typical PF fuel delivery 
system 1 may operate with an air to coal mass ratio of 1.5 to 2.5 with coal density of 1200 
to 1400 kg/m^ compared to 1 .2 kg/m^ for air at standard atmospheric conditions. Tlius, 
meeting the desired mass ratio results in a very dilute mixture of co4 on a yolumetric basis, 
on the order of one part in 1000 by volume. 

20 Assuming that the particles of coal are small enough and the acoustic jBrequencies 

and the firequencies of perturbations associated with the acoustics are low enough for the 
solid particles to exhibit negligible slip ( both steady and unsteady), the soimd speed can be 
assumed to be non-dispersive (that is constant with ftequency) and the volumetric phase 
fraction of the mixture could be determined through the Wood equation: 

25 

N 

1 

Pmix^lix f=l A^f 
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Including the effect of the compliance introduced by the conduit 12 (in this case a 
circular pipe of naodxjlus E, radius R and wall thickness t) 



1 1 . 
-..^ +cr where crs 



Pmix^ measured Pmix^mix 



Utilizing the relations above, the speed at which sound travels within the piping 
system of a representative coal / air mixtures is shown in Fig. 4 as a function of air /coal 

1 0 mass ratio. For liiis example, the pure air was assumed to have a density of 1 .2 kg/m^3 and 
a sound speed of 365.9 ru/s and the coal was assumed to have a density of 1400 kg/m'^S and 
a sound speed of 2439 m/s. As shown, the effect of increasing coal j&action, i.e. decreasing 
air/coal ratio is to decrease the sound speed. Physically, adding coal particles effectively 
mass loads the mixture, while not appreciably chaugmg the compressibility of the air. Over 

15 flie parameter range of interest, the relation between rnix^ sound speed and air / coal ratio 

is well behaved and monatonaic. 

While the calibration curves based on predictions from first principles are 
encouraging, using empirical data mapping from sound speed to air / coal ratio may result m 
improved accuracy of the present invention to measure the air/coal fractions of the mixture. 

20 However, it has been discovered that the physical properties of pulverized coal/air 

mixtures are generally such that there will be velocity slip at all but very low frequencies 
(on the order of <l-2 Hz for nominally 50 lam coal particles in air), as shown in Figs. 7 and 
8 which will be described in greater detail hereinafter. 

Fig. 5 shows the measured speed of sound as a function of frequency for an actual 

25 coal/air mixture 12. The sound speed was measured utilizing passive listening techniques 
of tbe present invention as described herein. The frequency dependence of the sound speed 
was determined by applying a Capon array-processiag algorithm at multiple narrow 
frequency ranges between 50-300 Hz thereby determining a frequency specific acoustic 
propagation velocity. In this particular example, the data was obtained wherein the coal/air 
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mixture was ftowing at nominally 100 ft/sec witii an air-to-coal mass ratio equal to 1 .8. The 
coal particles were nondnally 50 \xm in size, representative of pulverized coal typicaUy nsed 
in power generation and other industrial applications. A magnified view of the coal 
particles that were used for this test is shown in Fig. 3. 

Further shown in Fig. 5, the sound speed increases with increasing frequency and 
asymptotes toward a constaat value. The sound speed asymptote at higher frequencies is 
essentially the sound speed propagating through air only with no influence of the suspended 
particles. Also, it is apparent that the sound speed of the coal/air mixture has not reached the 
. quasi-steady limit at the lowest frequency for which sound speed was measured. The sound 
speed is continuing to decrease at the lower frequency limit An important discovery of the 
present invention is that the speed at which sound propagates through dUute particles 
suspended in a continuous fluid is said to be dispersive. As defined herein, the speed at 
which acoustic waves propagate through dispersive mixtures varies with frequency. 

Measuring the sound speed of a mixture 12 at progressively lower and lower 
fi-equencies becomes inherently less accurate as the total length of the array of pressure 
sensors 15 - 18 (Axaperaturc), wMch define Ihe aperature of the array, becomes small 
compared to the wavelength of the acoustics. In general, the aperture should be at least a 
significant fr action of a wavelength of the sound speed of interest In a particular 
embodiment soimd speed data was recorded with an array of four sensors, spaced at 12 
inches, for a total aperture of three feet At 50 Hz, a 1000 ft/sec sound wave has a 
wavelength of 20 ft. Thus, the aperture of this particular array (approx. 36 inches) spanned 
only 3/20flas of a wavelength, and the array's abihty to accurately resolve sound speeds 
below this was clearly impaired. It is an ioiportant aspect of the present invention that the 
abihty to resolve sound speed at low frequencies is directiy related to aperture of the array. 
Consequently longer arrays are used to resolve sound speeds at lower frequencies. As 
shown in Fig. 6, the standard deviation associated with determining the speed of sound in 
an is shown as a function of frequency for three arrays of varying aperture, namely 1,5 ft, 3 
fraud 10 ft. 

Given the practical constraints in accurately measuring sound speeds at ultra-low 
frequencies, the data suggests that utilizing a quasi-steady model to interpret the 
relationship between sound speed, measured at frequencies above those at which the quasi- 
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Steady model is applicable, and the air-to-fixel ratio would be problematic, and may, in fact, 
be impractical. Thus, the key to understanding and interpreting the composition of coal/air 
mixtures through sound speed measurements Ues in the dispersive characteristics of the 
coal/air mixture. 

5 In accordance with the present invention the dispersive nature of the system utilizes 

a fibrst principles model of the interaction between the air and particles. This model is 
viewed as being representative of a class of models that seek to account for dispersive 
effects. Other models cpuld be used to account for dispersive eflfects without altering the 
intent of this disclosure (for example, see the paper titled '"Viscous Attenuation of Acoustic 
1 0 Waves in Suspensions" by R.L. Gibson, Jr. and M.N. Toksoz), which is incorporated hereui 
by reference. The model allows for slip between the local velocity of the continuous fluid 
phase and that of the particles. The drag force on the particles by the continuous fluid is 
modeled by a force proportional to the difference between the local fluid velocity and that 
of the fluid particles and is balanced by inertial force: 

15 

where K - proportionality constant^ Uf = fluid velocity. Up - particle velocity, pp = particle 
density and Vp = particle volume. 

20 The effect of the force on the continuous fluid phase by the fluid particles is modeled as a 
force term in the axial momentum equation. The axial momentum equation for a control 
volume of area A and length Ax is given by: 



25 where P = pressure at locations x and Ax, ^ = volume ficaction of the particles, pf = fluid 
density. 

The particle drag force is given by: 
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where Cd = drag coefficient, Ap = frontal area of particle and pf = fluid density. 
5 Using Stokes law for drag on a sphere at low Reynold's number gives the drag coefBcient 



as: 



24 



TAjLt 



10 where Dp = particle diameter and |i = fluid viscosity. 



Solving for K in this model yields: 



15 



4 

Using the above relations and 1 -dimensional acoustic modeling techniques, the following 
relation can be derived for the dispersive behavior of an idealized fluid particle mixture. 



/ 



I 



1 + 



PpPp 



Pf 



1 + 0) 



2 2 



J 



20 



IQ the above relation, the fluid SOS, density (p) and viscosity (^) are those of the pure phase 
fluid. Vd is the volume of individual particles and <^p is the volumetric phase fraction of the 

particles in the mixture. 
25 Two parameters of primary interest in pulverized coal measurements are particle 

size and air^to-fuel mass ratio. To this end, it is of interest to examine the dispersive 
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characteristics of the mixture as a function of these two variables. Figs. 7 and 8 show the 
dispersive behavior for coal/air mixtures with parameters typical of those used in pulverized 
coal deliver systems. 

In particular Fig. 7 shows the predicted behavior for nominally 50 ^m size coal in 

5 air for a range of air-to-fuel ratios . As shown, the effect of air-to-fuel ratio is well defined 
in the low firequency limit. However, the eJBfect of the air-to-fuel ratio becomes 
indistinguishable at higher frequencies, approaching the sound speed of theptore air at high 
frequencies (above '-100 Hz). 

Similarly, Fig. 8 shows the predicted behavior for a coal/air mixture with an air-to- 

10 fuel ratio of 1.8 with varymg particle size. TMs figure illustrates thatparticle size has no 
influence on either the low frequency limit (quasi-steady ) sound speed, or on the high 
frequency limit of the sound speed. However, particle size does have a pronounced efifect 
in the transition region. • 

Figs. 7 and 8 iUustrate an important aspect of t^^ Namely, that 

15 the dispersive properties of dilute anixtures of particles suspended in a continuous fluid can 
be broadly classified into three frequency regimes: low JBrequency range, high frequency 
range and a ti*ansitional frequency range. As best shown in Fig. 8, the speed of sound 
propagatmg through the mixture is substantially the sa3aaie regardless of the particle size in 
the low frequency range. In the low frequency range the inixture exhibits a quasi-steady 

20 model or a no slip (non-dispersive) characteristic. As shown in the intermediate frequency 
range, the speed of sound propagating through the iiiixture is dependent on the size of the 
particle, and thus exhibits dispersive characteristics. For the high frequency range, the 
speed of sound propagating through the mixture is unaffected by the particles. In other 
words, the speed of sound in the higher frequency range propagating through the mixture is 

25 substantially equally to the speed of sound propagating through the just the fluid with the 
particles having for effect, which will be described in greater detail hereinafter. 

Knowmg the effect of dispersion on Ihe speed of sound through a mixture as 
described herein before, one will appreciate that to determine the concentration of the 
mixture (e.g., air/fuel ratio), the frequency of the measured acoustic wave is within the low 

30 frequency range that exhibits little or no slip (non-dispersive/quasi-steady state), as best 
shown in Fig. 7. Further, one will appreciate that to determine the particle size within the 
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mixture 12, the frequency of the measured acoustic wave is within the intermediate 
frequency range that exMbits dispersive characteristics, as shown in Fig. S. 

Although the effect of particle size and air-to-fuel ratio are inter-related, the 
predominant effect of air-to-fuel ratio is to determine the low frequency limit of the sound 
5 speed to be measured and the predommate effect of particle size is to determine the 
frequency range of the ti'ansitional regions. As particle size increases, the j&equency at 
which the dispersive properties appear decreases. For typical pulverized coal applications, 
this transitional region begins at fairly low frequencies, - 2Hz for 50 pm size particles. 

In the low frequency regime, f^^^ The 
10 j&equency range for which the no-slip, quasi-steady approximation is valid is a function of a 
variety of parameters including particle size, continuous phase viscosity, particle shape and 
particle density. 

The quasi-steady (no-slip condition) sound speed is given by the low firequency limit 
of the above relation, where AFR is air/fuel ratio: 



15 



1 



sa * 



AFR 



Note that particle size does not affect the low frequency limit of the sound speed. Referring 
to Fig. 9, the sound speed was measured using an embodiment of the present invention 

20 having eight sensors at 20.5 inch spacing, averaged from 20-40 Hz, for a range of air-to- 
coal mass ratios. The sound speed predicted for the coaFair mixtures using the quasi-steady 
model are also presented. As shown, although the general trend is captured, i.e. sound 
speed decreases with increased coal loading, the error is significant, rendering a first 
principle interpretation, based on a quasi-steady model inadequate. 

25 In the high firequency limit, the dispersion relation predicts the sound speed with 

■ 

asymptote towards the sound speed of the pure fluid. 



^mtxi^'=>'^)'=^ fluid 
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Interestingly, the higli frequency limit is independent of both particle size and air-to-fixel 
ratio. 

Given the difficulties measuring sufficiently low frequencies to apply the quasi- 
steady model and recognizing that the high frequency sound speed contains no direct 
iafoimation on either particle size or air-to-fuel ratio, it becomes apparent that the 
dispersive characteristics of the coal/air mixture should be utilized to determine particle size 
and air-to-ftiel ratio based on speed of sound measurements. 

As described hereinbefore, the flow meter 10 of tihe present invention includes the 
ability to accurately determine the average particle size of the coal in the PF/air nuxture 
within the pipe 14 and the air to fiiel ratio. Provided there is no appreciable slip between 
the air and the solid coal particle, the propagation of one dimensional soimd wave through 
multiphase mixtures is influenced by the eJBEective mass and the effective compressibility of 
the mixture. For an air transport system, the degree to which the no-slip assumption applies 
is a strong function of particle size and frequency. In the limit of small particles and low 
frequency, the no-slip assumption is vaUd. As the size of Hie particles increases and the 
frequency of the sound waves increase, the non-slip assumption becomes increasing less 
valid. For a given average coal particle size, the increase in slip with frequency causes 
dispersion, or, in other words, the somd speed of the naixture to change with frequency. 
With appropriate calibration the dispersive characteristic of a mixture will provide a 
measurement of the average particle size, as well as, the air to ftiel ratio (particle/fluid ratio) 
of the mixture. 

Using the model described above, which yields the equation shown below, and 
experimentally detennined sound speed as ftmction of frequency, the present invention 
includes an optiraization procedure to simultaneously determine particles size and AFR in 
coal / air mixtures: 



1 
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Referring to Fig. 10 there is shown an optimization procedure in accordance with the 
present invention in which the firee parameters of an analytical model are optiinized to 
ininimize an error fimction. For illustration purposes, the error function utOized is the sum 
of the differences of the sound speeds between an analytical model and the experimentally 
5 determinedsoundspeedas a function of frequency: ■ 

a; yX^(/)modei ^(/) measured)^ 
f-ftaw 

The results of the optimization procedure applied to data recorded from an array of sensors 

10 listening to flow in a six inch circular duct, 50 (xm particle size, 100 fit/sec air flow rate with 
anair-to-fuelratio of 1.8 is showninFig. 11. The measured and optimized--model- 
predicted soimd^eed is shown. As shown, the model captures the transitional frequency 
range well and provides a good estimate of the afr-rto-fiiel ratio. 

The results of Hie optimization procedure applied to a series of data sets with 

15 varying air-to-fiiel ratio is shown in Yig, 12. Note for this optimization the particle size was 
held constant over the range of data sets. 

As suggested herein before, the length of the array of pressure sensors should be at 
least a significant fraction of a wavelength of the sound speed of interest The significant 
fraction of the wavelength maybe at least 30 percent of the wavelength, however, this 

20 fraction may be less than 30 percent depending on the desired accuracy of the measurement, 
the measured wavelength and/or the strength of the acoustic wave (e.g., low signal/noise 
ratio). Therefore, the length of the array is dependent on the frequency of the sound speed 
of interest (fi-equency being inversely proportional to wavelength), wherein the frequency of 
the sound speed of interest is dependent on the measurement to be determined (e.g., 

25 air/particle ratio and particle size) and the dispersion characteristics of the mixture. For 
example, the low frequency range of the plot of the speed of sound (quasi-steady state) for 
measuring the concentration of the mixture (e.g., air/particle ratio), shown in Fig. 7, is lower 
as dispersion of the mixture increases. As described herein before, the dispersion 
characteristics of the mixture is dependent on the size of the particles among other factors. 

30 As the particle size increases, the dispersion becomes greater, and as the particle size 
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decreases, the dispersion becomes lower. Consequently, the length of the array is a function 
of the size of the particle within the nibctiire, and fherefo 
traxisition point (low frequency cut-off) between the low firequeiicy ra^^ 
ititermediate frequency decreases ia frequency as the particle size increases. 
5 For example when measuring the concentration of the mixture, as the size of the 

particles increase, the low frequency cut-off decreases and thus, the acoustic wavelength of 
interest increases to tihtereby necessitate the length of the array to be longer. Conversely, as 
the size of the particles decrease, the low frequency cut-off i^^ 

wavelength of interest decreases to thereby necessitate the length of the array to be shorter. 

1 0 Simply stated, the larger the particle, the longer the array and vice versa. The same 

comparison is true when determining the size of the particles within the mixture. However, 
for optimal performance of the flow meter, the measurement of tJie concentration of the 
mixture may require a longer array than the measurement of the particle size because 
measurement of the concentration is at a lower frequency (longer wavelength) than the 

15 intennediate fieqxiency (shorter wayelengfii) of the pa^ 

The lowest practical measurable frequency range is approximately 10-25 Hz, 
therefore the measurement of large particle may not be possible to measure the quasi-steady 
model, which may in some instances be less than 10 Hz (i.e., cut-off frequency less than 10 
Hz). Under these circumstances, the frequency of the speed of sound of interest is above the 

20 cut-off frequency. However, the measmred speed of sound is curve fit to a dispersion model 
of the inixture by varying the size of the particle and the composition of the mixture to 
determine the particle size and/or concentration of the mixture, as shown in Fig. 1 0 that will 
be described in greater detail hereinafter. 

While the length of the array is dependent on the particle size, the length may also 

25 be dependent on other parameters that define the amount of dispersion, such as mass of the 
particles and the viscosity of the fluid within the mixture. 

Another factor that defines (or effects) the length of the array of pressure sensors IS- 
IS includes the signal strength of the acoustic wave received by the processor. As the 
signal strength improves or is greater, the shorter the length of the array must be. The 

30 signal strength is dependent on a number of factors, such as the strength of the acoustic 
wave itself, the signal/noise ratio of the sensors, the matching of the sensors and others. 
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The spacing may be equi-spaced as shown in Fig. 1, however the flow meter 10 of 
the present invention contemplates that the sensors may have non-equal or uneven spacing 
therebetween. The sensors may be spaced any desired distance, provided the location or 
position of the sensors are known. For ported pressure sensors, the rninimum spacing is 
5 limited by mechaaical linaitatxons of the sensors. For strain-based sensors, such as PYDF 
bands described hereinafter, the compliance of the pipe limits the closeness of the spacings. 
For example, the more rigid the pipe, the greater the spacing of the sensors must be, and 
conversely, tilie more conapKant the pipe, the closer the sensors may be spaced. 

The spacing of the pressure sensors naay also be defined by the nurnber of sensors 

10 disposed within an array of a given length. The more sensors disposed within the array of a 
given length, the closer the spacing. The number of sensors disposed within an array is 
dependent on the required or desired accuracy of the flow meter 10, The greater the nimiber 
of sensors in the an*ay, a more precise meas\arement of the acoustic pressinre field can be 
achieved. Inother words, a greater number of samples or meas^ 

1 5 pressure wave over a given Imgth of the array (or wavelength) provided the sensors enable 
greater resolution iu the measurement of the acoustic wave to be measured or characterized. 

In addition to measuring the fluid to particle ratio of the naixture 12 and particle size 
within a pipe 14 using the measured speed of sound, the flow meter 1 0 further includes the 
ability to measure of volumetric flow rate of the mixture by comparing the difference of the 

20 speed of one dimensional sound waves propagating with and against the mean flow. 

This method of detemdning the volumetric flow rate of the particle/fluid mixture 12 
within pipe 14 relies on the interaction of the mean flow with the acoustic pressure field. 
The interaction results in sound waves propagating with the mean flow traveling at the 
speed of soimd (if the particle/liquid mixture were not flowing) plus the convection velocity 

25 and, conversely, sound waves traveling against the mean flow propagating at the speed of 
sound minus the convection velocity. That is, 

ClL^Clmix'' U 
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where aR = velocity of aright traveling acoustic wave relative to a stationary observer (i.e. 
the pipe 14), aL = velocity of a left trayeliag acoustic wave apparent to a stationary 
observer, amix fluid speed of sound (if the fluid were not flowing) and u == the mean flow 
velocity (assumed to be flowing flom left to right in this instance). Combiniag these two 
5 equations yields an equation for the meaa velocity. 

Therefore, by measuring the propagation velocity of acoustic waves in both directions 
.0 relative to the stationary pipe as described hereinbefore, the mean flow velocity can be 

calculated by multiplying the mean flow velocity by the cross-sectional area of the pipe 14. 

The practicality of using this method to determine the mean flow is predicated on 
tbe ability to resolve the sound speed in both directions with sufficient accuracy to 
determine the volumetric flow. For typical liquid measurements, flow velocities are 
15 typically at - 10 ft/sec and sound speeds of -4000 fl/sec. Thus axial mach numbers are on 
the order of 10/4000 of 0.0025: For a -17-10% accuracy in flow rate (+/- Ifl/sec), the sound 
speed of the upstream and downstream propagating waves would need to be resolved to +/- 

0.5/4000 or 1 part in 8,000. 

However, for PF/air mixture flows, axial flow velocities are nominally around 70 
20 ft/sec withno flow sound speeds of -700 ft/sec. This results in mach numbers of -0.1, 

approximately 2 orders of magnitude greater than typical liquid flows. For pulverized fiiel 
flows, to resolve the flow rate to 10% accuracy (or +/- 7 flt/sec), one would have to resolve 
the sound speed to +/- 3.5 ft/sec, or 3.5/700 or 1 part in 200. 

For the soimd speed measurement, the flow meter 10 utilizes simflar processing 
25 algorithms as those employed herein before. The temporal and spatial frequency content of 
sound propagating within the process piping 14 is related tlurough a dispersion relationship, 

a)=kanHx 

The wave number is Ic, which is defined as 'k^lid'k, o is the temporal firequency in rad/sec, 
and a„ux is the speed at which sound propagates within the process piping- For this cases 
30 where soimd propagates in both directions, the acoustic power is located along two acoustic 
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ridges, one for the soimd traveling with the flow at a speed of amix + Vmix and one for the 
sound traveling against the flow at a speed of amix - Vmix- 

The k-w plot shown in Fig. 13 illustrates the fimdamental principle behind sonar 
• based flow measure, namely that axial airay s of pressure sensors can be used in conjunction 

5 with sonar processing techniques to determine the speed at which naturally occurring 
turbulent eddies convect within a pipe. Fig. 13 shows a k-© plot generated for acoustic 
sound field of a coal/air mixture flowiag through a pipe. Two acoustic ridges are clearly 
evident. Each of the slopes of the two depicted acoustic ridges respectively defines the 
speed of sound traveling with and agamst the mean flow, respecti^^ A parametric 

10 optunization method was used to determine the "best" line representing the slope of the 
acoustic ridge. 

Further, Fig. 13 illustrates the abihty of the present invention to deteraiine the 
velocity of a fluid moving in a pipe. Fig. 14 shows a wavenumber-fiequency plot (k-w plot) 
of unsteady pressure. The contours represent the relative signal power at all combinations of 

15 frequency and wavenumber. The highest power "ridges" represent the acoustic wave with 
slope of the ridges equal to Ihe propagation speed. The dashed lines show the best-fit two- 
' variable maximization of the power with the two variables being sound speed andflow 
velocily. The right-side ridge represents the acoustic wave traveling in the same direction 
as the bulk flow and therefore its slope is steeper than the left-side ridge that represents the 

20 acoustic wave traveling in the opposite direction of the buUc flow. This indicates that the 
acoustic wave traveling in the same dhection of the flow is toveling faster than the acoustic 
wave traveling in the opposite direction of the bulk flow relative to the stationary sensors 
located on the pipe. 

The pressure sensors 15-18 described herein may be any type of pressure sensor, 
25 capable of measuring the unsteady (or ac or dynamic ) pressures within a pipe 14, such as 
piezoelectric, optical, capacitive, resistive (e.g., Wheatstone bridge), accelerometers (or 
geophones), velocity measuring devices, displacement measuring devices, etc. If optical 
pressure sensors are used, the sensors 15-18 may be Bragg grating based pressure sensors, 
such as that described in US Patent AppUcation, Serial No. 08/925,598, entitled " High 
30 Sensitivity Fiber Optic Pressure Sensor For Use In Harsh Environments", filed Sept. 8, 

1997, now U.S. Patent 6,016,702. Alternatively, the sensors 14 may be electrical or optical 
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Strain gages attached to or embedded in the outer or inner wall of the pipe which measure 
pipe waU strain, including microphones, hydrophones, or any other sensor capable of 
measuring the unsteady pressures within the pipe 14. In an embodiment of the present 
invention that utilizes fiber optics as the pressure sensors 14 they may be connected 
5 individually or may be multiplexed along one or more optical JSbers using wavelength 
division multiplexing (WDM), time division multiplexing (TDM), or any otiier optical 

multiplexing techniques. 

For any of the embodiments described herein, the pressure sensors, iricluding 
electrical strain gages, optical fibers and/or gratings among others as described herein, may 

10 be attached to the pipe by adhesive, glue, epoxy, tape or other suitable attachment means to 
ensure suitable contact between the sensor and the pipe 14. The sensors may alternatively 
be removable or permanentiy attached via Imown mechanical tech^ 
mechanical fastener, spring loaded, clamped, claua shell arrangement, strapping or other 
equivalents. Alternatively, the strain gages, including optical fibers and/or gratings, may be 

15 embedded in a composite pipe. If desired, for certain applications, the gratings may be 
detached firom (or strain or acoustically isolated firom) the pipe 14 if desired. 

It is also w:ithin the scope of the present invention th^^ 
technique may be used to measure the variations ia strain ia the pipe, such as highly 
sensitive piezoelectric, electronic or electric, strain gages attached to or embedded in the 

20 pipe 14. 

In certain embodiments of the present invpntion^ a piezo-electronic pressure 
transducer may be used as one or more of the pressure sensors 15-18 and it may measure 
the unsteady (or dynamic or ac) pressure variations inside the pipe 14 by measuring the 
pressure levels inside of the pipe. In an embodiment of the present invention, the sensors 14 

25 comprise pressure sensors manufactured by PCB Piezotronics. In one pressure sensor there 
are integrated circuit piezoelectric voltage mode-type sensors that feature built-in 
microelectronic amplifiers, and convert the high-impedance charge into a low-impedance 
voltage output SpecificaUy, a Model 106B manufactured by PCB Piezotronics is used 
which is a high sensitivity, acceleration compensated integrated circuit piezoelectric quartz 

30 pressure sensor suitable for measuring low pressure acoustic phenomena in hydrauUc and 
pneumatic systems. It has the unique capabiUty to measure small pressure changes of less 



wo 03/091671 



25 



PCT/US03/12956 



than 0.001 psi imder high static conditions. The 106B has a 300 mV/psi seBsitivity and a 
resolution of 91 dB (0.0001 psi). 

The pressixre sensors incorporate a built-in MOSFET microelectronic amplifier to 
convert the high-impedance charge output into a low-impedance voltage signal. The sensor 

5 is powered &om a constant-current source and can operate over long coaxial or ribbon cable 
without signal degradation. The low-impedance voltage signd is not affected by 
triboelectric cable noise or insulation resistance-degrading contaminants. Power to operate 
integrated circuit piezoelectric sensors generally takes the form of a low-cost, 24 to 27 
VDC, 2 to 20 mA constant-current supply. A data acquisition system of the present 

10 invention may incorporate constant-cuirent power for directly powermg integrated cii'cmt 

piezoelectric sensors. 

Most piezoelectric pressure sensors are constructed with either compression mode 
quartz crystals preloaded ia a rigid housiag, or unconstraiued tourmaline crystals. These 
designs give the sensors microsecond response times and resonant frequencies in the 
15 hundreds of kHz, with minhnal overshoot or ringing. Small diaphragm diameters ensure 

spatial resolution of narrow Shockwaves. / ; 

The output characteristic of piezoelectric pressure sensor systems is that of an AC- 
coupled system, where repetitive signals decay imtil there is an equal area above and below 
the origmal base line. As magnitude leyels of the monitored event fluctuate, the output 
20 remains stabilized around the base line with the positive and negative areas of the curve 
remaining equal. 

Furthermore the present iavention contemplates that each of the pressure sensors IS- 
IS of the flow meters 10,70 may include a piezoelectric sensor 104 -107 that provides a 
piezoelectric material 110 to measure the unsteady pressures of the fluid/particle mixture 12 
25 as shown in Fig. 20. The piezoelectric material, such as the polymer, polled 

fluoropolymer, polyvinylidene fluoride (PVDF), measures the strain induced within the 
process pipe 14 due to unsteady pressure variations within the process mixture 12. Strain 
within the pipe is transduced to an output voltage or current by the attached piezoelectric 
sensors 104-107. 

30 As best shown in Fig. 21, the PVDF material 1 10 is adhered to the outer surface of a 

steel strap 112 that extends around and clamps onto the outer surface of the pipe 14. The 
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piezoelectric sensing element is typically confonnal to allow complete or nearly complete 
circtimferential measurement of induced strain. The sensors can be formed jBrom P VDF 
films, co-polymer jBilms, or flexible PZT sensors, similar to that described in 'Tiezo Fikn 
Sensors technical Manual" provided by Measurement Specialties, Inc., which is 
incorporated herein by reference. The advantages of this technique are the following: 

1. Non«intrusive flow rate measurements 

2. Low cost 

3 . Measurement technique requires no excitation source. Ambient flow noise is 
used as a source. 

4. Flejcible piezoelectric sensors can be mounted in a variety of configurations to 
enhance signal detection schemes. These configurations ijaclude a) co-located 
sensors, b) segmented sensors with opposing polarily configurations, c) w 

sensors to enhance acoustic signal detection a^^ 

detection, d) tailored sensor geometries to miniinize sensitiyity to pipe modes, e) 
differencing of sensors to eUminate acoustic noise fiom vortical signals. 
5- Higher Temperatures (140C) (co-polymers) 

While the present invention illustrates that the array of pressure sensors comprises a 
plurality of like sensors, the present invention contemplates that any combmation pf 
different or similar pressure sensors may be used wiflm 

While tibie present invention is capable of measuring solid particles suspended in a 
fluid, one wiU appreciate that other multi-phase mixture^ or flows may be measured using 
an array of sensors, such as steam flow. It is further recognize the that effects of dispersion 
on large solid particles hi a fluid would be sunilar to large droplets of a liquid dispersed in a 
gas or air, and thus similar considerations when measuring the steam quality and droplet 

size should be addressed. 

It should be understood that any of the features, characteristics, alternatives or 
modifications described regarding a particular embodiment herein may also be applied, 
used, or incorporated with any other embodunent described herein. 

Although the invention has been described and illustrated witib respect to exemplary 
embodiments thereof, the foregoing and various other additions and omissions may be made 
therein and thereto without departing from the spirit and scope of the present invention. 
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Claims 

What is claimed is: 

' 1 . An apparatus for measuring at least one parameter of a particle/fluid mixture in a 

5 pipe, comprising: 

a spatial array of at least two pressure sensors, disposed at different axial locations 
along the pipe, and each measnring an unsteady pressure within the pipe at a corresponding 
axial location, each of said sensors providing a pressure signal indicative of the unsteady 
pressure within the pipe at said axial location of a corresponding one of said sensors; and 
10 a signal processor, respopsive to said pressiire signals, w^ 

indicative of the at least one parameter of the mixture in the pipe. 

2. The apparatus of claim 1 wherein each sensor measures an acoustic pressure and 
provides a signal indicative of an acoustic noise within the pipe. 

15 

3. The apparatus of claim 1, wherein the signal processor, responsive to said pressure 
signals, provides a signal indicative of the speed of sound propagating through the mixture 
in the pipe, . 

20 4. The apparatus of claim 3 wherein said signal processor comprises logic which 
calculates a speed at which sound propagates along said spatial array. 

5. The apparatus of claim 3 wherein said signal processor comprises logic which 
calculates a frequency based signal for each of said acoustic pressure signals. 

25 

6. The apparatus of claim 4 wherein said acoustic pressure signals each comprise a 
frequency based signal and wherein said signal processor comprises logic which calculates 
a ratio of two of said frequency based signals. 

30 7. The apparatus of claim 1 comprising at least three of said sensors. 



wo 03/091671 



28 



PCT/US03/12956 



8. The apparatus of claim 3 wherein the signal processor comprises logic which 
calculates a fliiid composition of the rnixture in the pipe. 

9. The apparatus of claim 1 wherein at least one of said pressure sensors measures a 
5 circumference-averaged pressiure at said axial location of said sensor. 

1 0. The apparatus of claim 9 wherein at least one of said pressure sensors includes a 
piezoelectric sheet material. 

10 11. The apparatus of claim 1 wherein the piezoelectric sheet material is polarized 
fluoropolymer, polyvinyhdene fluoride (PVDF). 

12. The apparatus of claim 1 wherein at least one of said pressmre sensors measures 
strain on the pipe. 

15 

13. The apparatus of claim 3 wherein the frequency based sound speed is determined 
utilizing a dispersion model to determine the at least one parameter of the inixture. 

14. The apparatus of claim 3 whereia the array of acoustic sensors are spaced 

20 sufficiently such that the entire length of the array is at least a significant fraction of the 
measured wavelength of the acoustic waves being measured. 

15. A method for measuring at least one parameter of a particle/fluid mixture in a pipe, 

said method comprising: 
25 measuring unsteady pressures within the pipe at at least two predetermined axial 

measurement locations along the pipe to provide a pressure signal indicative of the unsteady 
pressure within the pipe at each of the at least two predetermiaed axial measurement 
locations; and 

calculating the at least one parameter of the particle/fluid mixture in the pipe using 
30 the unsteady pressiire measured at the axial measurement locations. 
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25. The method of claim 24 wherein the measured unsteady pressures are acoustic 
pressures to provide a signal indicative of an acoustic noise within the pipe. 

26. The method of claim 25, wherein the calculating the at least one parameter uses an 
5 acoustic pressure to calculate a speed of soimd propagating in the pipe. 

27. The method of claim 26 wherein the calculating the at least one parameter uses an 
acoustic pressure to calculate a speed at which sound propagates along said spatial array. 

10 28. The method of claim 26 wherein the calculating the at least one parameter uses an 
acoustic pressure to calculate a frequency based signal for each of said acoustic pressure 
signals. 

29. The method of clakn 27 whereiQ said acoustic pressure signals each comprise a 

1 5 frequency based signal and whereia said signal processor comprises logic which calculates 
a ratio of two of said frequency based signals. 

30. The method of claim 24 wherein measuring the unsteady pressures is at at least three 
of said sensors. 

20 

3 1 . The method of claim 26 wherein the calculating the at least one parameter uses an 
acoustic pressure to calculate a fluid composition of the mixture in the pipe. 

32. The method of claim 24 wherein measuring unsteady pressure iacludes measuring a 
25 circimiference-averaged pressure at at least an axial location of a sensor. 

33. The method of claim 32 wherein measuring unsteady pressures uses at least one of 
said pressure sensors includes a piezoelectric sheet material. 

3 0 34. The method of claim 24 wherein the piezoelectric sheet material is polarized 
fluoropolymer, poly vinylidene fluoride (PVDF). 



wo 03/091671 PCT/US03/12956 

30 



34. The method of claim 24 whereia at least one of said pressure sensors measures strain 
on the pipe. 

5 35. The method of claim 26 wherein the frequency based sound speed is determined 
utilizing a dispersion model to deterniine the at least one parameter of the im 



10 



36. The method of claim 26 wherein the array of acoustic sensors are spaced sufficiently 
such that the entire length of the array is at le^t a significant fraction of the measured 
wavelength of the acoustic waves being measured. 
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